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The ultracold neutron (UCN) upscattering rate from para- 
deuterium contaminated solid D2 near T — K is calcu- 
lated and found to be proportional to the atomic fraction 
of para contaminant (x par a) as t^ p ss 700 x para sec -1 . With 
Xpara — 0.33, corresponding to room temperature equilib- 
rium, the upscatter time is 4.4 msec. The implication is 
that the para-deuterium contaminant is the temperature- 
independent energy source responsible for the observed up- 
scattering rate in the Los Alamos spallation driven UCN 
source. 



I. INTRODUCTION 

The possibility to use low neutron absorption cold ma- 
terials, such as solid deuterium, as a moderated source 
of ultracold neutrons (UCN) has been long discussed [ 1 1 . 
However, as pointed out by Golub and Pendlebury [2|, 
one can in fact attain a much higher UCN density than 
implied by the thermal equilibrium arguments given in 
0; if we consider UCN produced by inelastic down- 
scattering (in energy) of cold or thermal neutrons, the 
reverse process can be relatively suppressed. In the case 
of moderator with two energy levels, a ground state and 
an excited state separated by an energy A, the upscat- 
ter rate of a UCN with energy E ucn and the downscatter 
rates of a neutron with energy E ucn + A rates are related 
according to microscopic reversibility: 



Odo 



-A/kT 



(1) 



The implication of this result is that, at very low temper- 
atures, UCN can be accumulated in a material until the 
rate of production P equals the rate of loss, 1/t, due to 
nuclear absorption or leakage from a storage volume that 
surrounds the converter material; the density in equilib- 
rium is 



p = Pt 



(2) 



This is the so-called "superthermal" UCN production 
process, first described by Golub and Pendlebury in rela- 
tion to thier proposed superfluid 4 He UCN source. 00] 
The theory of the superthermal process for superfluid ^rle 
has been now very well established, both by observation 



of upscattered UCN from a storage volume pj, and by 
direct measurement of the UCN density in a nearly losses 
superfluid 4 He filled magnetic trap. 

The possibility of the use of other materials of low neu- 
tron absorption cross section as a superthermal source 
was first investigated by Golub et al. || after Mr. Kiliv- 
ington observed the production of UCN in a prototype 
superfluid 4 He source, operated at the Institut Laue- 
Langevin, even when it was empty of He; the expla- 
nation by Golub was that there was a film of dirt on 
the inner surface of the UCN converter/storage vessel 
and this dirt was operating as a superthermal source. 
This led to the proposal of the "thin film superthermal 
source" : Unlike the case of 4 He, all other materials have 
some neutron absorption cross section. For materials like 
solid D2, the absorption lifetime is 120 ms, compared to 
the effective lifetime of about 880 sec of a UCN in su- 
perfluid 4 He at low temperature, in which case /3-decay 
is the ultimate limiting loss mechanism. The idea of the 
thin film source is that the UCN density in a storage vol- 
ume of volume V that contains some converter material 
is independent of the volume of the converter material V c 
because the total production rate is proportional to V c , 
while the net effective lifetime r is proportional to 1/V C 
(in the limit where upscattering and absorption in the 
converter material dominates all other loss mechanisms), 
and V c cancels in Eq. 0. V c determines the time constant 
for filling the volume V to an equilibrium density. 

The figure of merit for a superthermal source is given 
by Pt for the converter material. Superfluid 4 He is at 
least an order of magnitude better than any other mate- 
rial; in this case, even though P is relatively low, the long 
intrinsic lifetime more than makes up the difference in 
production rate compared to other materials. However, 
there are situations where use of a less efficient material 
is warranted. For example, a thin film source might be 
better suited near the core of a nuclear reactor because 
less converter material is required to achieve the intrinsic 
maximum UCN density. 

Another situation where a material with a large P but 
short t was first pointed out by Pokotilovski J7|] in con- 
nection with pulsed neutron sources, such as spallation 
sources. The idea is that the converter material is con- 
nected to a UCN storage volume only for a short time 
during and after the neutron pulse, after which a valve 
is closed, separating the conversion and storage regions. 
In this case, one can enjoy a relatively large P together 



with an net effectively large r. Serebrov elaborated upon 
this idea with the proposal of a "UCN Factory" ; follow- 
ing this lead, we developed a prototype UCN source that 
turns out to follow the original Pokotilovski idea more 
closely, but ours incorporates a cold neutron flux trap, 
and vertical UCN extraction from the solid D2 converter 
to counteract the solid deuterium positive potential step. 
In our recent studies of solid deuterium as a spalla- 
tion driven generator of ultracold neutrons, the apparent 
production rate is roughly an order of magnitude less 
than expected, and there is a lack of temperature de- 
pendence on the net production rate, contrary to the 
prediction by Golub et al. ||[| The reason for the re- 
duction appears primarily due to a loss of UCN before 
they can exit the solid deuterium. By elimination of var- 
ious loss sources within the material that can account 
for the temperature independent loss, we came to the 
conclusion that the UCN loss is likely due to upscatter- 
ing from para-deuterium contaminant with solid (ortho-) 
deuterium converter material; in addition, anecdotal ev- 
idence for this was offered by Serebrov. The following 
calculation estimates the upscatter rate and was moti- 
vation for incorporating a para/ ortho converter [pj . The 
preliminary results indicate the following calculation is 
correct; the data and analysis will be presented in a full 
publication by the collaboration. We felt it important, 
however, to make available our initial estimate of the up- 
scatter rate due to the para contaminant. 



under normal conditions, and we can assume that the 
room temperature composition persists as the D2 sample 
is cooled and solidified, and is rather permanent in our 
experiment. 

We can assume the sample is at absolute zero; we 
therefore can neglect upscattering from the vibrations of 
the lattice which represents a different upscatter mech- 
anism, roughly independent of upscattering from the 
para contaminant. Furthermore, we can make an ap- 
proximation that the molecules are fixed in the lattice 
and neglect recoil effects; this approximation is valid be- 
cause a) Eq < On, the Debye temperature, b) the mo- 
mentum transfer is in the form of angular momentum 
from the molecular rotation to the outgoing wave, c) 
the UCN energy Ei « Od, making direct excitations 
of the lattice unlikely. Finally, we assume that all the 
para molecules are in the J = 1 rotational state, while 
the ortho molecules arc in the J — state. This final 
assumption is valid because the rotational levels, unlike 
the nuclear spin state, are closely coupled to the thermal 
bath, as mentioned above. 

Normally, we assume that P— wave scattering is sup- 
pressed for cold and ultracold neutrons. However, in this 
case, the scattering is from a molecule; the separation 
between the deuterium atoms in the molecule is compa- 
rable to the wavelength of an outgoing upscattered UCN 
(with final energy 7 meV) and we do not have the usual 
P— wave scattering suppression found in nuclear scatter- 
ing from a single isolated nucleus. 



II. PROPERTIES OF D 2 AND ASSUMPTIONS IN 
THE CALCULATION 



III. UPSCATTER RATE 



The deuterium molecule exists in the ground state as a 
total nuclear spin S = and rotational angular momen- 
tum J = 0, and is referred to as ortho— deuterium. 

The first excited state has 5 = 1 and J = 1, and the 
energy relative to the ground state is 



E = 7 meV. 



(3) 



We refer to this state as para— deuterium. 

Of course, at high temperature, the S, J states are oc- 
cupied according to their degeneracy and Boltzmann fac- 
tor; just considering the nuclear states, the ortho-states 
(S = 0, 2) and para-states (S — 1) are occupied at high 
temperatures according to their spin degeneracies: 



1 + 3 + 5 



[J. GO] %ortho 



1 + 5 
1 + 3 + 5 



0.67 



(4) 



because the degeneracy is given by 2S* + 1. 

The rotational states are coupled strongly to the lattice 
and we can assume the J populations follow the Boltz- 
mann distribution with rapid equilibration. On the other 
hand, the transitions between ortho and para is very slow 



Much of Young and Koppel (YK) || is irrelevant to 
this calculation and in addition is not entirely transpar- 
ent. Our problem is simple enough to begin from first, 
very basic, principles. We will see than only a single ma- 
trix element from YK is required for our calculation, as 
subject to the above constraints and approximations. 

Starting with Fermi's Golden Rule for transitions to a 
continuous state (see fLOfl, Sec. 43) 



dw }l = {2ir/h)\U fi \ 2 5(E f - EJdvf 



(5) 



where dvjfi is the transition rate between the levels /, i, 
Avf is the density of final states, and Ufi is the matrix 
element of the Hamiltonian between the initial and final 

states, 



U(r) 



2irk 2 



(axSir- Fi) + a 2 S(r- r 2 )\ 



(6) 



where 1, 2 label the two atoms of the molecule, m is the 
neutron mass, and a\ 2 are the scattering lengths, chosen 
to represent the ortho- or para- states as required. 

First, let us take the incoming and outgoing neutron 
states as plane waves (as did YK). The incident (UCN) 



wave must represent the UCN density p u in a normaliza- 
tion volume V. Thus, 



1 



Vv = —f=\/ 9uV e %kr = sfptyi 



(7) 



where the factor 1/yV normalizes the incident wave to 
the volume V, the second factor gives the total UCN 
n = p u V within V, the subscript i,n means incident 
wave, normalized, and ipi conforms to the states used by 
YK. 

The outgoing wave can be written 



The density of final states kf is 

Vdkf _ Vkjdkfdn 



dv 



I 



(2tt) s 



(2tt) 3 



(8) 



(9) 



In Fermi's Golden rule, the (5-function ensures energy 
conservation. In the present case, 



E f -Ei = h 2 kj/2m - h 2 k 2 /2m - E 



(10) 



and we have <5(/(fc/)). Recall the well-known relation- 
ship, 



*(/(*)) 



1 



l/'(*o)| 



S(x-x ) 



where xq is the zero of f(x). A further assumption fcj = 
Ei = is valid for UCN in this case. Then the momentum 
transfer k — kf, in YK's notation. 
Therefore, the upscatter rate is 



dwfi 



2tt /2nfif 



h \ m 
x |(^ / ,J , = 0,S" = l|^;a# i J = l I 5 = 0)|' 

VkUkt 



h% 



x 8(kf - yf2mE /n 2 )-^- 

which must be integrated over kf, this is easily done, and 
we write the matrix element | ^ Oi| 2 , yielding a total rate, 
taking Vx para pu 2 scattering molecules in V, 



h sr-^ 

Wfi,tot = ^Vx para pD 2 Pucn — k \ } a t 



where 



and 



k a = J2mE h 2 = 1.83 x 10 8 cm _1 



-k = vq = 1.2 X 10 5 cm/s. 



(11) 



(12) 



(13) 



We further note that Wfijot/V — p ucn so we have a loss 
rate of 



T = iTTX para p D2 V Q \^ a i\ 



(14) 



Thankfully, YK have calculated the matrix elements. 
We only need 

\(k f , J' = 0, S' = 1| J^ ai \k t = 0, J = 1, S = 0)| 2 = (15) 

?> ' ■' ' (16) 

(17) 



^^(1011:00)14,11' 

T a inc\Aoi\ 



from Eq. (A5) of KY and Appendix B. Also, (a = .742A 
so fc a/2 = 0.68) 

\A 2 m \ = \2i 3l (k Q a/2)\ 2 = 4 x j^koa/2) 2 = 4 x (.21) 2 - .18 

(18) 

(19) 

(20) 



so 



2 
inc 



\J2 a *\ 2 = °- i3a 

Taking a inc — 4.04 x 10~ 13 cm gives 

r = 3.2 X 1Q- 20 XparaPD^ 1 



and, with x para = 0.33 corresponding to room tempera- 
ture equilibrium, and po 2 =3x 10 22 /cm 3 gives 



T = 316 sec -1 ; r = 3.2 



(21) 



Finally, zero-point motion of the lattice leads to a 
smearing-out of the final state, and a reduction in the 
cross section. This effect is given by the Debye- Waller 
factor, which for solid deuterium is a factor of 0.72 for 
a momentum transfer corresponding to 7 meV. |11| The 
net upscatter rate and lifetime is thus 

r = 221sec~ 1 ; r = 4.4 msec. (22) 



IV. DISCUSSION 

The important implication of this calculation is that 
the upscatter rate from ortho-deuterium is independent 
of temperature and persists at absolute zero. The agree- 
ment between this simplified calculation and experiment 
is satisfactory and appears to have identified a principal 
limitation to the present Los Alamos UCN source density. 

The surprising result is that the P— wave scattering is 
significant; the outgoing plane wave state carries away 
the rotational angular momentum of the molecule. This 
is possible because the impact parameter, corresponding 
to the molecular size, is comparable to the outgoing neu- 
tron wavelength, as is illustrated in Eq. |18L 



V. CONCLUSION 

Contamination of a solid deuterium UCN source by 
the para-molecular state will severely limit the attainable 
UCN density. Given the stability of the para-deuterium 
state (0.02%/hr conversion rate in the solid) it is clear 
that the gas must be converted before freezing. The de- 
tails of the Los Alamos converter will be published, but it 
follows standard techniques. |q] Our preliminary results 
support the above calculated rate. 

Our result has serious implication for solid deuterium 
sources operated close to the core of a nuclear reactor. 
In this situation, the gamma flux will be sufficient to 
drive molecular excitations, and the ortho-para ratio will 
likely be near the high temperature equilibrium. This 
implies that solid deuterium might be useless in such an 
environment, although the presence of free charge in the 
lattice might catalyze the conversion. We are presently 
investigating these and other effects expected in a high 
radiation environment. We further point out that the 
radiation levels in our proposed spallation driven source 
are many orders of magnitude less than a reactor so we 
expect no problems in this regard. 

A more complete calculation of the upscattering rate 
and related issues is in progress by the authors. 
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